Core hole-electron correlation in coherently coupled molecules 
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We study the core hole-electron correlation in coherently coupled molecules by energy dispersive 
near edge X-ray absorption fine-structure spectroscopy. In a transient phase, which exists during 
the transition between two bulk arrangements, 1,4,5,8-naphthalene-tetracarboxylicacid-dianhydride 
multilayer films exhibit peculiar changes of the line shape and energy position of the X-ray ab- 
sorption signal at the C _R"-edge with respect to the bulk and gas phase spectra. By a comparison 
to a theoretical model based on a coupling of transition dipoles, which is established for optical 
absorption, we demonstrate that the observed spectroscopic differences can be explained by an in- 
termolecular delocalized core hole-electron pair. By applying this model we can furthermore quantify 
the coherence length of the delocalized core-exciton. 



X-ray absorption measurements at high-resolution 
beam lines of third-generation synchrotron sources are 
a powerful tool for the investigation of the electronic and 
vibronic structure of molecules [IH3]- From an analysis 
of high-resolution data one can gain insight into the com- 
plex response of the electronic system and of the molec- 
ular frame on the electronic excitation. Upon excitation 
of a core electron into an unoccupied molecular orbital 
a particularly strong reaction of the electronic system in 
the vicinity of the core excitation occurs. This core hole- 
electron correlation, which is often referred to as core 
exciton, has to be included for a correct description of 
the experimental data. However, the nature of core exci- 
tons is not clear in molecular materials, since experimen- 
tal information, e.g. on the spatial expansion, is usually 
complicated to extract from X-ray absorption data due 
to various other interfering effects (4j-[6] . Thus it is not 
clear whether the excitation is mainly localized on one 
molecule or on a subgroup of a molecule, or if it extends 
over several molecules. 

In this work we demonstrate the importance of the 
core hole-electron correlation in near edge X-ray absorp- 
tion fine structure (NEXAFS) data of organic molecules. 
On the example of the organic dye 1,4,5,8-naphthalene- 
tetracarboxylicacid-dianhydride (NTCDA) we show that 
under certain circumstances the molecules exhibit an ex- 
traordinarily strong intermolecular coupling, which leads 
to a red shift of electronic transitions and a substan- 
tial narrowing of the vibronic progressions in the NEX- 
AFS spectra with respect to the gas phase. The respec- 
tive experimental data can be described by a coupling 
of the transition dipoles to the neighbouring molecules 
thus resulting in a significant intermolecular delocalisa- 
tion of the core hole-electron pair. While this effect is 
well-established for optical excitations |S] , the calcu- 
lated transition dipoles in core-excitations are substan- 
tially smaller and thus a more local character of the ex- 
citation may generally be expected. 

The delocalisation of the core-excition upon X-ray ab- 



sorption implies not only the transfer of the excited 
electron but also of the core hole. While for diatomic 
molecules, such as iV 2 and O2, experimental evidence for 
the non-local character of the core hole exists [9TfTT]. an 
intermolecular delocalisation of a core hole is not easily 
understood. 

For the model system NTCDA (see molecular struc- 
ture in Fig. [lb) on a Ag(lll) surface different struc- 
tural phases exist for multilayer films [51 [T2J [TB] • Growth 
at substrate temperatures below 100 K results in mor- 
phologically smooth films with preferentially flat lying 
molecules and without long range order [12j . Heating 
the samples to temperatures above 180 K leads to an in- 
creasing crystallinity in the films and changes the molec- 
ular orientation to a more upright standing configuration. 
During this transition a transient phase occurs, which ex- 
ists for only about 10 min and which shows the character- 
istic changes of the signature of the NEXAFS spectra of 
particular interest for this work. The fast acquisition of 
X-ray absorption data is thus crucial and the experiments 
were therefore performed at the beamline UE52-PGM of 
BESSY II in Berlin. The beamline was operated in the 
energy dispersive mode, where the beam is focused on 
the sample without an exit slit, thus leading to a spatial 
dispersion of photon energies on the sample in vertical 
direction [2]. By use of an electron spectrometer with 
adapted spatial resolution (VG-SCIENTA R4000) NEX- 
AFS spectra can be recorded in a multi-channeling mode 
with a possible acquisition time of less than one second 
and an energy resolution better than 50 meV at the C PC- 
edge. 

r igure [TJl) shows the C K-NEXAFS spectra of a four 
monolayer (4 ML) thick NTCDA film which was prepared 
and measured at a substrate temperature of 95 K. The 
absorption signal of the it* -resonances between 282 eV 
and 290 eV shows a characteristic dependence on the po- 
larization of the incident light, which allows determining 
the average molecular orientation angle with respect to 
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Figure 1. a) C K-NEXAFS spectra of a NTCDA multilayer 
film as prepared at 95 K, recorded with s- (black line) and 
p-polarization (grey line) of the incident beam, b) Close 
up of the 7r*-resonances (denoted as A-C) recorded with p- 
polarisation, compared to the result of a Hartree-Fock cal- 
culation with the GSCF3 code [15] . The theoretical energy 
scale was shifted and scaled to fit the experimental data. The 
calculated electronic transitions are marked according to the 
respective carbon ls-sites CI to C4 (see inlay) and final state 
orbitals n* , 2-n* , 3-7T*. The shaded region indicates the en- 
ergy window recorded in the energy dispersive mode, which 
is displayed in panel c). 



the substrate surface [16]. The strong dichroism in Fig- 
ure [IJi) indicates a preferential flat lying orientation of 
the NTCDA molecules in the film, consistent with pre- 
vious investigations [31 [T5]. Figure [IJj) shows the -n*- 
resonances of the p-polarized spectrum on an expanded 
energy scale. According to Hartree-Fock calculations [T5] 
plotted on the bottom of Figure [Tja) the spectrum can 
be explained by the superposition of several electronic 
transitions from the four symmetry-inequivalent carbon 
Is sites CI to C4 (see inset in Figure [l}))) to the final 
state orbitals l7r*, 2tt* , 3tt* . In addition, the electronic 
transitions show a characteristic coupling to vibrations, 
leading to the rich fine structure of the three main sig- 
nals, denoted by A, B, and C in Figuref!]}) [3J. While the 
spectra in Figure [iji) and b) were recorded by scanning 
the photon energy in the normal beamline mode with an 
acquisition time of about 30 min per spectrum, Figure 
[TJ;) shows an energy-dispersive NEXAFS experiment of 
the energy regime of feature B (shaded region in Figure 
[lb)). The integration time was 30 seconds and the data 



is in very good agreement with feature B of Figure [lp) . 
The fine structure is well resolved and according to the 
calculations in Figure [TJd) the signal in the photon en- 
ergy range between 285.1 and 286.0 eV corresponds to 
the two electronic transitions C4 ls-27r* and C3 ls-27r*, 
which are originating from carbon atoms located on the 
naphthalene core (see Figure [T|d)). As demonstrated in 
Ref. [3J these electronic transitions lead to the excita- 
tion of vibrations of the aromatic core with vibrational 
energies of around 90 meV. 

The colour plot in Figure [2^i) shows the evolution of 
feature B with time during an energy dispersive NEX- 
AFS experiment. The sample (thickness about 4 ML) 
was prepared at 95 K and then quickly (in about 6 min) 
heated to 220 K in order to induce the transition from 
flat lying to upright standing molecules. On the y-axis 
t=0 refers to the time when a constant sample tempera- 
ture of 220 K was reached. On the bottom of Fig. [2^,), 
i.e. directly after preparation at 95 K, the spectra re- 
semble the typical absorption signal for films prepared at 
low temperature (see Fig. [TJ;)). This is demonstrated by 
Fig. [2]d) , which displays spectra derived after integration 
of the intensity map in Fig. [2^) over intervals of 3 min 
at the times indicated by the white horizontal lines. 

Shortly after a substrate temperature of 220 K is 
reached the line shape of feature B is significantly 
altered. Starting at t—2 min the intensity of the 
absorption peaks at 285.4 eV and 285.7 eV increases, 
while it decreases in-between. In addition, the line width 
of the most prominent signals decreases, which can be 
observed best in the spectrum in Fig. J2Jd) extracted 
after 8 min at 220 K. Moreover, an energy shift of the 
main spectral components towards lower energy occurs, 
which is indicated by the white dashed lines in Fig. 
[2]}). The energy shift sums up to about ~ 80 meV after 
35 min While the overall intensity increases until about 
t=12 min, the absorption signal is constantly decreasing 
after the sample was kept at 220 K for more than 12 min. 

This variation of general intensity can be explained by 
the change of orientation angle of the NTCDA molecules 
from fiat lying to a more upright configuration, which oc- 
curs after the sample temperature is increased to 220 K 
[12) . This leads to an immediate intensity increase in 
the C K-NEAXFS data recorded with s-polarisation [TB] . 
However, this intensity increase is counteracted by a 
change of film morphology, since the increased sample 
temperature also leads to a roughening of the originally 
closed film and to formation of 3D islands. Due to the 
limited probing depth of about 1 nm [TB] of the Auger 
electrons mainly contributing to the NEXAFS signal in 
electron yield mode, this second process leads to an over- 
all reduction of the absorption signal. Moreover, the mor- 
phology change includes material transport and occurs 
on a slower time scale than the change of orientation an- 
gle, thus accounting for the observed delayed intensity 
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a) Orientation of the dipole moment for 
coupled molecules 
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Figure 2. a): 2D colour map of the energy dispersive NEX- 
AFS spectra of feature B recorded with s-polarized light plot- 
ted against time. Yellow refers to high, blue to low intensity. 
The 4 ML thick film was prepared at 95 K and heated to 
220 K within 1 min. t=0 refers to the time when a constant 
sample temperature of 220 K was reached. The vertical white 
dashed lines are guidelines for the eye to illustrate the energy 
shift of the main components, b) NEXAFS spectra of feature 
B derived by integration of the intensity map of Figure a) over 
3 min at the times indicated by the horizontal white dashed 
lines (black dots: original data; solid line: 3-point average). 



decrease [19] . 

We will in the following concentrate on the partic- 
ular changes of the spectral signature of the NTCDA 
C K-NEXAFS which occur during the transition. Fig- 
ure [3]i) compares the NEXAFS spectra of feature B for 
four different NTCDA samples. On the bottom gas phase 
data is displayed, which is representative for the isolated 
molecule [17] . Film growth at low sample temperature 
( T=95 K) results in the disordered phase with preferen- 
tially flat lying molecules (second spectrum from bottom) 
[3], while annealing to temperatures above 180 K results 
in the ordered phase with more upright molecular ori- 
entation (third spectrum from bottom) [3] . The data of 
the transient phase during the transition is plotted on 
top. The line shape in the disordered phase resembles 
the gas phase, while the entire spectrum is red-shifted by 
about 100 meV upon condensation. The ordering of the 
NTCDA molecules upon annealing leads to a slight de- 
crease of the line width and to a more pronounced shoul- 
der on the low energy side of feature B. This trend is 
followed if the transient phase is examined. The most 
intense peak at 285.7 eV is now very sharp and the low 
energy shoulder has developed into a well distinguishable 
peak at 285.4 eV. Moreover, additional peaks can be re- 
solved at 285.5 eV and 285.6 eV as well as at 285.8 eV. 
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Figure 3. a) Illustration of the coupling of NTCDA molecules 
in the transient phase. The transition dipole moments are 
oriented perpendicular to the molecular plane, b) Schematic 
energy level diagram of adjacent NTCDA molecules after res- 
onant core-excitation. See text for details, c) Result of a 
coherent exciton scattering calculation: The adiabatic com- 
ponent of the electronic excitation (0-0) and the correspond- 
ing vibronic progression (0-1, 0-2, 0-3) are plotted for isolated 
(gray) and coupled molecules (black), d) NEXAFS spectra of 
NTCDA feature B in the gas-phase |17| . in the disordered 
condensed phase !3i, in the ordered condensed phase after an- 
nealing [3], and in the transient phase (from bottom). The 
spectra were normalized arbitrarily for better comparison of 
the line shape. 



The energetic separation of these signals from the promi- 
nent peaks at 285.4 eV and 285.7 eV corresponds well 
to the vibronic energy of the NTCDA modes coupled to 
the C4 ls-27r* and C3 ls-27r*- transitions, respectively 
[3J. We thus straightforwardly interpret the spectrum of 
feature B in the transient phase as due to the prominent 
adiabatic transitions C4 ls-27r* at 285.4 eV and C3 ls- 
2-7T* at 285.7 eV with adjacent vibronic states. The main 
difference to the other condensed and gas phase spectra 
lies in a reduced line width of the involved signals and in 
a change of the vibronic envelope, which obviously shifts 
intensity towards the adiabatic transition in the transient 
phase. 

An explanation for this observation can be given on 
the basis of a coupling of the transition dipole moments 
of adjacent molecules. The coupling strength strongly 
depends on the arrangement of the molecules and on the 
symmetry of the respective transition moments. In the 
dipole approximation the coupling strength V nm between 
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molecule n and m can be described by 
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were /2„ and pi m are the transition dipole moments on 
molecule n and m, respectively, which are separated by 
f. For a C ls-27r* transition with a dipole moment per- 
pendicular to the molecular plane the coupling is thus 
largest for a configuration as sketched in Figure [3^). 

The proposed dipole-dipole coupling mechanism is il- 
lustrated in Figure [3]}). The C ls-2-7r* excitation of 
molecule A couples to the same transition on the neigh- 
boring molecule B. 

The changes of the vibrational envelope observed in 
our NEXAFS data can be understood in the framework 
of the coherent exciton scattering (CES) theory [20T - 
123"] . This one-particle Green's function approach obeys 
Dyson's equation and describes the excitation of coher- 
ently delocalized quasi-particles. In other words this 
means that the hole-electron pair propagates to the ad- 
jacent molecule fast compared to the molecular vibration 
frequencies. Therefore the originally excited molecule re- 
mains preferentially in its vibrational ground state. As 
a result, intensity is shifted towards the adiabatic transi- 
tion in the vibrational envelope observed in experiment. 
This is illustrated by Figure pp) , which shows the result 
of a CES calculation of the vibronic progression of a sin- 
gle molecule and of coupled molecules in a configuration 
as sketched in Figure j3ji) 24 . The intensity ratio of 
the adiabatic transition (0-0) to the higher vibrational 
states (0-1, 0-2, and 0-3) changes strongly, and in case of 
sufficiently strong coupling the oscillator strength accu- 
mulates in the 0-0 transition. Moreover, a marked energy 
shift of the entire vibronic progression towards lower en- 
ergy can be observed in case of the coupled molecules, 
which also matches our experimental observation quali- 
tatively (see Figure [2^)). 

In the transient phase we also obverse a reduced line 
width of the C ls-27r* transitions, as evidently seen in 
Figure [3]i) . For a quantitative analysis the full width 
at half maximum (FWHM) of the gas-phase and of the 
transient phase spectra was estimated to 211 (±10) meV 
and 74 (±10) meV respectively, for the most prominent 
signals in the data of [3]i) . According to Knapp et al. [23J 
the decrease of the FWHM A of the absorption signal 
with the number of coherently coupled molecules N can 
be described by the equation 
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With the FWHM derived from our NEXAFS data of 
isolated and coherently coupled molecules, i.e. Aj SO ;. 



and A 



coup . ) 



respectively, a number of about N—8±3 can 



be estimated in the transient phase. 



In conclusion, we have reported on a X-ray absorption 
experiment which demonstrates the importance of a cor- 
rect consideration of the core hole- electron correlation 
for condensates of organic molecules. During the transi- 
tion from a disordered to an ordered film arrangement, 
a transient phase occurs for NTCDA/Ag(lll) which is 
characterised by a distinctly different spectral signature. 
The spectra of this transient phase exhibit a reduced 
line width, a red shift, as well as a substantial change 
of the vibronic profile. All these experimental observa- 
tions can be explained by the coupling of the transition 
dipole moments of adjacent molecules, a model, which 
is well-established for optical spectroscopies. This re- 
sult has two immediate implications: Firstly, it shows 
that the excitation delocalises over several molecules (a 
number of about 10 can be determined from our data). 
This comprises the existence of a dclocalised core hole, 
which is involved in the excitation. This implication is 
not understood and needs further consideration. Note 
that other models have been considered but could not 
explain our experimental observation. 



Secondly, the very accurate explanation of our exper- 
imental data by this model based on the coupling of 
transition dipoles strongly suggests that the oscillator 
strength of the transition dipole moments involved in 
the core-to-valence excitations are substantially larger 
than estimated on the basis of calculations on isolated 
molecules. This discrepancy may be explained by the 
additional contribution of the wave function overlap be- 
tween neighboring molecules [26-- 29 . For close-packed 
7r-conjugated molecules in a face-to-face arrangement 
neighboring dyes can have a considerable overlap of the 
frontier orbital wave functions [30--33J. This can lead to 
increased intermolecular charge transfer transition dipole 
moments [M] [35] ■ This so-called short-range component 
of the coupling becomes significant, e.g., for naphthalene 
dimers already at a separation of about 6 A|36j and can 
thus be expected to be substantial in the present case, 
were the intermolecular distance between the NTCDA 
molecules is only around 3 A. Further efforts to develop 
a comprehensive theoretical description of the core exci- 
tation in molecular condensates are thus highly desirable. 
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